Excitation of MHD waves in magnetized anisotropic cosmologies 
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The excitation of cosmological perturbations in an anisotropic cosmological model and in the pres- 
ence of a homogeneous magnetic field was studied, using the resistive magnetohydrodynamic (MHD) 
equations. We have shown that fast-magnetosonic modes, propagating normal to the magnetic field 
grow exponentially and saturated at high values, due to the resistivity. We also demonstrate that 
the jeans-like instabilities are enhanced inside a resistive and the formation of condensations formed 
within an anisotropic fluid influence the growing magnetosonic waves. 

PACS numbers: 



I. INTRODUCTION 

Magnetic fields are known to have a widespread pres- 
ence in our Universe, being a common property of the 
intergalactic medium in galaxy clusters (Kronberg 1994), 
while, reports on Faraday rotation imply significant mag- 
netic fields in condensations at high redshifts (Kronberg, 
Perry & Zukowski 1992). Studies of large-scale magnetic 
fields and their potential implications for the formation 
and the evolution of the observed structures, have been 
the subject of continuous theoretical investigation (sec 
Thorne 1967, Jacobs 1968, Ruzmaikina & Ruzmaikin 
1971, Wasserman 1978, Zel'dovich, Ruzmaikin & SokolofT 
1983, Adams, Danielsson & Rubinstein 1996, Barrow. 
Fcrrcira & Silk 1997, Tsagas & Barrow 1997, Jedamzik, 
Katalinic & Olinto 1998, Barrow, Maartens & Tsagas 
2006, for a representative, though incomplete list). Mag- 
netic fields observed in galaxies and galaxy clusters are 
in energy equipartition with the gas and the cosmic rays 
(Wolfe, Lanzetta & Oren 1992). The origin of these 
fields, which can be astrophysical, cosmological or both, 
remains an unresolved issue. 

If magnetism has a cosmological origin, as observations 
of /^G fields in galaxy clusters and high-redshift proto- 
galaxics seem to suggest, it could have affected the evolu- 
tion of the Universe (Giovannini 2004, Barrow. Maartens 
& Tsagas 2006). There are several scenarios for the gen- 
eration of primordial magnetic fields (e.g. see Grasso 
& Rubinstein 2001). Most of the early treatments were 
Newtonian, with the relativistic studies making a rela- 
tively recent appearance in the literature. A common 
factor between almost all the approaches is the use of the 
MHD approximation, namely the assumption that the 
magnetic field is frozen into an effectively infinitely con- 
ductive cosmic medium (i.e. a fiuid of zero resistivity). 
With a few exceptions (e.g. sec Fennelly 1980, Jedamzik, 
Katalinic & Olinto 2000, Vlahos, Tsagas & Papadopou- 
los 2005), the role of kinetic viscosity and the possibility 
of non-zero resistivity have been ignored. Nevertheless, 
these aspects are essential for putting together a compre- 
hensive picture of the magnetic behavior, particularly as 
regards the non-linear regime. The electric fields associ- 



ated with the resistivity can become the source for par- 
ticle acceleration, while the induced non-linear currents 
may react back on the magnetic field (Vlahos, Tsagas & 
Papadopoulos 2005). 

Many recent studies have used a Newtonian or a Fried- 
mann - Robertson - Walker (FRW) model to represent 
the evolving Universe and super-imposed a large-scale 
ordered magnetic field. The magnetic field is assumed 
to be too weak to destroy the FRW isotropy and the 
anisotropy, induced by it, is treated as a perturbation 
(Ruzmaikina & Ruzmaikin 1971, Tsagas & Barrow 1997, 
Durrer, Kahniashvil & Yates 1998). Current observa- 
tions give a strong motivation for the adoption of a FRW 
model but the uncertainties on the cosmological Stan- 
dard Model are several. Therefore, the limits of the ap- 
proximations and the effects one may lose by neglecting 
the anisotropy of the background magnetic field, should 
be investigated. Within this context, the formation of 
small-scale structures and the excitation of resistive in- 
stabilities in Bianchi-Type models has been explored sev- 
eral years ago (Fennelly 1980). Nevertheless, the excita- 
tion of MHD-waves in curved spacetime and their sub- 
sequent temporal evolution, is far from being exhausted 
(Papadopoulos, Vlahos & Esposito 2001). 

In the present article we explore the evolution of a 
magnetized resistive plasma in an anisotropic cosmolog- 
ical model. We begin with a uniform plasma driving the 
dynamics of the curved spacetime (the so-called zeroth- 
order solution). This dynamical system is subsequently 
perturbed by small-scale fluctuations and we study their 
interaction with the anisotropic background, searching 
for imprints on the temporal evolution of the perturba- 
tions' amplitude. In particular: 

In Section II, we present the system of the field equa- 
tions appropriate to describe the model under considera- 
tion. In Section III, we solve this system analytically, to 
derive the zeroth-order solution. Accordingly, in Section 
IV, we extract the first-order perturbed equations. In 
Section V, we derive the dispersion relation for the mag- 
netized cosmological perturbations and in Section VI, we 
perform a numerical study of their evolution, using a 
fifth-order Runge - Kutta - Fehlberg temporal integration 
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scheme. Finally, in Section VII, a perturbation analysis 
over purely gravitational fluctuations reveals an inherent 
Jeans-like instability. 

Our results suggest that in a resistive plasma, within 
an interval of lO^^sec after the beginning of the interac- 
tion process, fast-magnetosonic modes are excited, grow- 
ing exponentially in time and is saturated at high values. 
In this way, magnetic field perturbations can be retained 
at large amplitudes, for sufficiently long time-intervals 
{At ~ 10^^ sec), resulting in the enhancement of the am- 
bient magnetic field (dynamo effect). In addition, the 
resistive plasma enhance the condensations that can be 
formed within the anisotropic fluid due to a gravitational 
instability, which, in turn, influence the growth of the 
magnetosonic waves. 



II. THE FIELD EQUATIONS 

We consider an axisymmetric Bianchi-Type I cosmo- 
logical model, driven by an anisotropic and resistive per- 
fect fluid, in the presence of a time-dependent magnetic 
field, B = B{t)x. The corresponding line-element is writ- 
ten in the form 

ds^ = -df + R^{t)dx^ + S'^{t)[dy'^ + dz^] (1) 

The evolution of a curved spacetime in the presence of 
matter and an e/m field, is determined by the gravita- 
tional field equations 



7^ 



1 



(2) 



(in the system of units where h = 1 = c) , together with 
the energy-momentum conservation law 



T^'^ = 
and Maxwell's equations 



(3) 

(4) 
(5) 



In Eqs (2) - (5), Greek indices refer to the four- 
dimensional spacetime (in accordance, Latin indices refer 
to the three-dimensional spatial section) and the semi- 
colon denotes covariant derivative. Furthermore, 7?.^^ 
and TZ are the Ricci tensor and the scalar curvature with 
respect to the background metric g^i,, while G is New- 
ton's gravitational constant. Eventually, F'"' is the an- 
tisymmetric tensor of the e/m field and is the corre- 
sponding current density. 

The energy-momentum tensor involved, consists of two 
parts; namely 



^ ~ ■'fluid "1" -'em 



(6) 



The first part, is due to an anisotropic perfect fluid 
source, of the form 



(7) 



where, p{t) is the energy density, Pi{t) are the compo- 
nents of the anisotropic pressure and the axial symme- 
try of the metric (1) implies that P2{t) = Psit)- Fi- 
nally, = dx^/ds is the fluid's four- velocity, satisfy- 
ing the conditions u^u^ = 
h^'^ = 9^" + u^u'^ being the projection tensor 
The second part, is due to the ambient e/m field 



1 and h^'^'u^ = 0, with 



em 



1 

47r 



1 



{F^'^F'^^g^p _ g^^'^F^pF^^) 



(8) 



where, for E ~ Q and B / / x, the non-zero components 
of the Faraday tensor in the curved spacetime (1) read 
(see Appendix A) 



pa; 



(9) 



As regards the current density J^, it may be determined 
by the invariant form of Ohm's law 



(10) 



where, Ug is the locally measured charge density and 77 
is the (finite) electric resistivity, in units of time. As a 
consequence of the Maxwell equations, we obtain J.*^ = 0. 
Assuming that that the fluid has zero net-charge, i.e. 
Ue = 0, Eq (10) reduces to J'^ = ^F^'^Uy. A vanishing 
net charge indicates that the perfect fluid consists, at 
least, of two components. 



III. THE BACKGROUND SOLUTION 

We look for axisymmetric Bianchi-Type I cosmologi- 
cal solutions to the Einstein-Maxwell equations (Appen- 
dices A and B), representing the background metric of 
our problem. In this case, Eqs (2) reduce to 




8TTGp{t) + GB'^{t) 



8nGpi{t) -GB^{t) 



(11) 





-1- 


(-1 






\RS 



(the dot denotes time-derivative) and Eqs (4), (5) yield 

dt[S^B{t)] = (12) 

Eq (12) has a clear physical interpretation: The magnetic 
flux through a comoving surface normal to the direction 
of the magnetic field, is conserved. 

On the other hand, the continuity equation (3) results 
in (Appendix C) 

W0 + f^] + |[..(t)-£] + 2fM0 + g] 







(13) 
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and the particles' number conservation law, reads 



IV. THE COSMOLOGICAL PERTURBATIONS 



(14) 
The system of Eqs (11) - (14) admits the exact solution 



to IQ 

Pl(t) ^'lo(v)^ P2(0 = -P2o(v)^ 



(15) 



t 



t 



where, the index "0" stands for the corresponding values 
at t = to a-iid to marks the beginning of the interaction 
between magnetized plasma and curved spacetime. Solu- 
tion (15), represents an anisotropic cosmological model, 
in which, the large-scale anisotropy along the x-axis, is 
due to the presence of an ambient magnetic field. The 
combination of Eqs (11) and (15) indicates that, initially, 
the total energy density is given by 



Po 



Bl 



(16) 



and the difference between fluid's pressure and the pres- 
sure of the magnetic field along the two anisotropic spa- 
tial directions is equal 



For any dynamical system, much can be learnt by in- 
vestigating the possible modes of small-amplitude oscilla- 
tions or waves. A plasma is physically much more com- 
plicated than an ideal gas, especially when there is an 
externally applied magnetic field. As a result, a variety 
of small-scale perturbations may appear. We first as- 
sume a uniform magnetized plasma in curved spacetime 
as background, which is perturbed by small scale fluctu- 
ations. In this article, the evolution of the background is 
described by the solution (15). 

Accordingly, we introduce first-order perturbations to 
the Einstein-Maxwell equations, by decomposing the 
physical variables of the fluid as 



p{t,z) = p{t) +Sp{t,z) 

Px{t,z) = pi{t) 

Py{t,z) = P2{t) - Sp{t,z) 

Pz{t,z) = P2{t) +Sp{t,z) 



(21) 
(22) 



and we insert the perturbed values (21) and (22) into Eqs 
(11) - (14), neglecting all terms higher or equal than the 
second order. The pressure perturbation Sp{t, z) intro- 
duces a longitudinal acoustic mode, propagating along 
the i-direction and therefore 



Sp{t, z) = Cl 5p{t, z) 



(23) 



P2a 



Stt 

Bl 

Stt 



1 



(17) 



where, Cg is the speed of sound. The four-velocity of 
the plasma fluid is perturbed around its comoving value, 

u^' = (1,0,0,0), as 

uf'it, z) = {l + (Su°(t, z), 0, 0, Su^{t, z)) (24) 



Eqs (17) lead us to identify 



Pw = Po P20 



(18) 



i.e. initially, when i?(to) = S(to), the two components 
of the anisotropic pressure were equal in absolute value, 
something that is confirmed also by Eq (13). Further- 
more, with the aid of Eqs (16) and (17), we obtain 



3 B"^ 
Pio + 2p2o = 7(po + 6-^) 
5 Stt 



(19) 



which, according to Eq (18), results in the equation of 
state for the matter-energy content at t = to; namely. 



Po 



lipo 



4^ 



(20) 



For i?o = 0, i.e. as regards the perfect fluid itself, we 
obtain that, initially, po = \po- Since po < our 
model corresponds to a semi-realistic cosmological model 
of Bianchi Type I. These models are crude, first order 
approximations to the actual Universe when we use cur- 
rently available theories and observations (Jacobs 1969). 



Then, the condition 
der, implies 



= —1, to the first leading or- 



5u°(t, z) = 



and, therefore, u^(t, z) 



(5w^(t,z). 



(25) 

Accordingly, 



p(t, z)u^{t, z) = p{t)5u^{t, z) + O2. 

As regards the perturbations of the e/m field, we con- 
sider that they correspond to a transverse e/m wave, 
propagating along the z-axis (kj j z)\ namely, 

E{t,z) = 8Ey{t,z)y (26) 
i3(t,z) = B{t)x + SB''{t,z)x (27) 

Therefore, now, the non-zero components of the Faraday 
tensor in curved spacetime arc modified as follows 



F"2 = i<5£;J'(t,z) = -f2o 

F^' = ^[B{t)+6B'\t,z)] = -F^' 



(28) 



In what follows, we take into account the so-called 
Cowling approximation (Cowling 1941), admitting that 
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59iiy = 0. Therefore, the evolution of the perturbed 
quantities is governed only by the energy-momentum ten- 
sor conservation, together with Maxwell's equations. 

To begin with, we perturb the particles' number con- 
servation law: Accordingly, Eq (C3), yields 

d,{5p) + pmXSu^ + ^p{jl + 2§) = ^ (29) 

We continue with Maxwell's equations: Then, from Eq 
(B2), using Eqs (21), (22), (24) and (25), we obtain 



ATT-\5Ey + SB{t)5u''] 

■n 



(30) 



Now, Eq (B3) becomes 

dt{S^5B^) - Sd,{5Ey) = (31) 
Eventually, the conservation equation (C2) results in 

dt[p{t)5u''- 



while, to the first leading order, Eq (CI) collapses to an 
identity. Eqs (29) - (32) arc the linearly independent 
first order perturbed Einstein-Maxwell equations in the 
curved background (1), which we intend to discuss. No- 
tice that, in the flat spacctime - zero resistivity limit, they 
reduce to Eqs (10.53a), (10.9) and (10.53c) of (Jackson 
1975), respectively. 

To develop the theory of small-amplitude waves in 
curved spacetime, we search for solutions to the linearized 
equations (29) - (32), in which all perturbation quantities 
are proportional to the exponential 



exp [i{kz 



ujdt)] 



(33) 



following the so-called adiabatic approximation 
(Zcl'dovich 1979, Birrell & Davies 1982, Padman- 
abhan 1993). In this context, the (slowly varying) 
time-dependent frequency of the wave is defined by the 
eikonal 



through the relation 



n = 



jdt 



(34) 



(35) 



Notice that, in Eq (33), z is the comoving coordinate 
along the i-axis and k is the comoving wave-number. In 
an expanding Universe, the corresponding physical quan- 
tities are defined as Zpfi = zS{t) and fcp/j ~ k/S{t), so 
that kpfiZpfi = kz. 



V. THE DISPERSION RELATION 

Before discussing the temporal evolution of the pertur- 
bation quantities, it is important to trace what kind of 
waveforms are admitted by this system, in the first place. 
In order to do so, we have to derive their dispersion rela- 
tion, D(k,Lu) = 0, at i = tQ. Provided that certain kinds 
of modes (such as acoustic, magnetosonic etc) do exist 
in the first place, they can be excited through their in- 
teraction with the anisotropic spacetime. An additional 
excitation, due to the non-zero resistivity, is also possible 
(Fennelly 1980). 

Accordingly, we assume a wave-like expansion for the 
perturbation quantities, of the form 

6p = A^e^C^^-/* Su' = Aue'^'''-^''^'^*'^ {36) 

= yi^e*^''^'-^''^'^*', = ylBe'(''^--/'*"''*)(37) 

Sp = yipe*(fc^-r-'rft) ^ c2^^g*(fc^-r'^<it) (38) 

Although the background quantities depend on time, in 
the search for a dispersion relation at t = to, we treat 
the perturbation amplitudes (Ats) as constants. In this 
way, our search for potential waveforms at t = to, is not 
disturbed by the inherent non-linearity introduced for 
t > to. Nevertheless, once the potential waveforms are 
determined, their interaction with the curved spacetime 
in the presence of an external magnetic field, implies that, 
for t > to, the time-dependence of their amplitudes is a 
priori expected. Using Eqs (36), Eq (29) is written in 
the form 

{Hr + 2Hs - iuj)5p ~ikp(t)6u'' (39) 
where, we have set 



R S 
Hft = — and Hs ~ — 
" R ^ S 

Furthermore, using Eqs (37), Eq (30) reduce to 



(40) 



Att k Att 

[iio - {Hr + Hs + —)]SEy = ~i-6B^ + —SB(t)Su'' 
rj b 7] 

(41) 

while, Eq (31) becomes 



{2Hs - iio)6B'' = i^5Ey 

J 



(42) 



Finally, Eq (32) yields 



p{t) + {Hr + 2Hs - iuj)p{t)]5u'' + i^ClSp = 



i 1 |,i|l,.,H„.. 

- i—kB{t)SB'' 
Att 



(43) 



With the aid of Eqs (12) and (14), the combination of 
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Eqs (39) - (43) results in 



which yields 



[ -u^ + kl^Cl - iuj{Hi, + 2Hs)] X 

4:7T 

= —u\{HR + 2Hs-iuj) X 

[ (HR-Hs-iu;){2Hs-iu;) + klf,] (44) 

where, u\ = (Bq/Attpq) is the (dimcnsionlcss) Alfvcn 
velocity. Eq (44) is the dispersion relation which deter- 
mines the possible waveforms admitted by this dynamical 
system for every t>to- 

We have to point out that, co, as defined by Eqs (34) 
and (35), has the usual meaning of the angular frequency 
of an oscillating process only in the short-wavelength 
(high-frequency) regime of the mode k (Mukhanov, Feld- 
man & Brandenberger 1992). In other words, the wave 
description in curved spacetime makes sense only when 
the physical wavelength along the direction of propaga- 
tion [Xph = AS'(t)] is much smaller than the correspond- 
ing horizon length [ins = ^5^(^)1; 



(52) 



Hs 



(45) 



Eq (45) implies that, in the anisotropic background (1), 
the wave description makes sense as long as 



uj ,kph~> Hr, Hs 



(46) 



for every t > Iq. In this limit, Eq (44) becomes surpris- 
ingly transparent, namely 

{LO^-kl,,C!)iLO^-kl,,)+zu;^ [c.2(l + ul) - kl^iCl + u\)\ 

. . 

Vanishing of the real part results in acoustic (a; = kp^Cs) 
and e/m (lj = kjjh) waves, while, vanishing of the imagi- 
nary part results in fast-magnetosonic waves 



This result has a clear physical interpretation: Every 
very-high-frequency perturbation of the dynamical sys- 
tem under study, is suppressed due to the finite resistiv- 
ity. Therefore, the only modes that survive in a resis- 
tive cosmological model, are the (low-frequency) MHD 
modes. In the next Section, we intend to discuss the 
evolution of these modes. 



VI. NUMERICAL STUDY OF THE MHD MODE 

In order to study the temporal evolution of the magne- 
tosonic modes for t > to , we assume that their amplitudes 
are no longer time-independent 

SEy = i^(t)e*('=^-^*'^^*), ,55^ = i?(t)e'('=^--/'*'^'^*(54) 
5p = p(t)e^('=^-/''^'^*) = C2p(t)e'('="--/'''^''*) (55) 

In Eqs (53) - (55), the wave- number k is related to the 
frequency lo through Eq (48) and, once again, we have 
taken into account the equation of state for the perfect 
fluid. 

We decompose the time-dependent amplitude of the 
perturbations (53) - (55) into a real and an imaginary 
part, as 



~p{t) = PR{t)+ipi{t) 

u{t) = UR{t) + iui{t) 

E{t) = ER{t)+iEi{t) 

B{t) = BR{t) + iBi{t) 



(56) 



(57) 



u:\\ + u\) = kUCl+u\) 



(48) 



something that reduces Eqs (29) - (32) to the following 
first order system 



In fact, in the zero-resistivity limit (ideal plasma), the 
obvious modes expected are the magnetosonic modes, 
which we recover. On the other hand, in most astro- 
physical situations we have (Jackson 1975) 



u\ « Cl 
In this case, Eq (47) reads 

.47r 



(49) 



{Lu^-kl,,Cl){oj^ + i—oj-kl,,) = Q (50) 



PR + ujpi - p^{ — fkui + -pR = 

t 2 
pi - LopR + po{ — )^kuR + -pi = 

Er + LuEi + kC-^f'^Bi + %Er + Att-Er 
t Zt rj 

+4^is„(^)i/2^^^0 
?7 t 



Ei - ojEr - kC-jf'^BR + ^Ei + Att-Ei + 



3 

2i' 



According to Eq (50), in the very high frequency limit 
where no acoustic waves are admitted, we are left with a 
waveform governed by the dispersion relation 



2 47r 2 
U) + i — UJ — kpf^ — 



(51) 



+ATr-BoC-^y/^ui = 
1] t 

Br + ujBi + jBr + kCjf/^Ei = 
Bi - luBr + ^Bj - kC-^Y'^ER = 



(58) 
(59) 

(60) 

(61) 
(62) 
(63) 
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1 



1, 



AtTIq 1] t 

-BluR~kCl{^)pi=Q 
to 



1 



PqUi - LUpoUR + (— — + 



47rto '? t 



--B^o^i + kCl{^)pR = 



(64) 



(65) 



Wc integrate numerically the system (58) - (65), using a 
fifth order Runge-Kutta-Fehleberg scheme with variable 
integration step. The time is measured in units of and, 
therefore, r = ^ > 1. In terms of r, the physical wave- 
number reads kph = k/y/r and the Hubble parameter 
along the yz-plane is written in the form Hs = {2t^)~^ . 
According to Eq (45), for a certain value of r, a wave is 
well inside the horizon as long as 



k > 



1 



(66) 



The validity of Eq (66) for long r- intervals determines the 
appropriate values of the comoving wave-number. Now, 
the analysis depends on where do we place the initial 
time, to- 

According to the Standard Model (Kolb & Turner 
1990), after nucleosynthesis, the Universe goes on ex- 
panding and cooling but nothing of great interest takes 
place until t ^ l[y^^ sec. At that time, the temperature 
drops to the point where electrons and nuclei can form 
stable atoms (recombination). Before that time, dur- 
ing the so-called radiation epochs photons couple strongly 
with matter, the main constituent of which is in the form 
of plasma. Therefore, the latest time at which plasma 
could play a role of cosmological significance is the recom- 
bination time {tR = 1.2 X 10^^ sec). In the limiting case 
where to = tR, the condition (66) reads fc ^ x 10~^* 

sec~^ and, therefore, an appropriate choice for k would 
be fc = 10-12 sec-^. 

In order to decide on the initial values of the unper- 
turbed quantities, we write Eq (16) in ordinary units, 
namely 



he' 



(67) 



We adopt a typical behavior for the energy-density, valid 
at the late stages of the radiation epoch (see Weinberg 
1972, Eq 15.6.42) 



poc^ = 1.45aT^ 



(68) 



where, T is the temperature and a is the black-body 
constant. At the time of recombination {to = 1.2 x lO^"^ 
sec, T = 4000 °K), we obtain pq = pqc'^ = 2.8 erg/cm^, 
which, through Eq (67), is effectively a choice on Bq; 
namely i?o — 7 gauss. Notice that this value lies barely 
within limits of the constraint 



2 Bq 



(69) 



a necessary condition to retain the anisotropy of the met- 
ric (Thorne 1967). Extrapolation of this result, along the 
lines of Eq (15), to the present epoch {tp ~ 15 x 10^ y), 
suggests that, today, the corresponding magnetic field 
should be Bp ~ 6.6 x 10""'^'' gauss. This value lies 
within limits of the upper bound for the present-day 
magnetic field strength, arising from the large-angular 
scale anisotropy of the microwave radiation background 
(MRB) at last scattering (Barrow, Ferreira & Silk 1997, 
Barrow, Maartens & Tsagas 2006) 



i? < 4 X 10 gauss 



(70) 



In fact, now, we may proceed even further, to estimate 
the amount of distortion which the expansion anisotropy 
along the x-axis (caused by the unperturbed magnetic 
field) induces to the microwave pattern at the present 
epoch. The contribution of a large-scale coherent mag- 
netic field to the microwave quadrupole anisotropy is 
given by (Madsen 1989) 



AT 



f>2 

o7re„ 



(71) 



where. Bp and ep denote the present values of the mag- 
netic field and the background radiation energy-density, 
while, z is the redshift at which the anisotropy begins 
to grow (in our case, at the recombination time, where 
z ~ 1100). The present value of the microwave back- 
ground temperature is Tp = 2.8 °K, corresponding to 
an energy-density of ep ~ 4.7 x 10""'^^ erg/cm^ for the 
radiation field. Accordingly, our analysis suggests that 
the present-day quadrupole anisotropy along the x-axis 
should be 



AT 



4.06 X 10" 



(72) 



i.e. almost four times larger than the corresponding 
COBE result. 

Taking into account that, initially, the unperturbed 
quantities are of the order of unity, wc normalize all the 
perturbation quantities at t = to, to 0.01 in cgs units. 
On the other hand, initially, the equation of state for the 
perfect fluid admits = 0.2, while, as regards the resis- 
tivity, we adopt the Spitzer relation (Krall & Trivielpiece 
1973) 



?/ = 10 



-2 (2L)3/^ 



sec 



(73) 



In a radiation-dominated background, we have (Kolb & 
Turner 1990) 



T 10^ 



(74) 



and therefore, during recombination, Eq (64) results in 
rj — 0.0645 sec. In order to demonstrate how rj may 
trigger instabilities, we consider three cases, namely: rj = 
0.0645 sec, ry = 0.0745 sec and tj = 0.0870 sec. 
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The output of the numerical integration consists of the 
electric and the magnetic field perturbations' amplitude 



0.035 



\SEy\ 

\SB'=\ 



El 



(75) 
(76) 



and illustrates their temporal evolution. In Fig. 1, we 
present the magnetic field perturbation versus time. We 
consider two cases: 

• For rj = (ideal plasma), the magnetic field per- 
turbation grows steeply at early times. It appears 
that, the interaction of the perturbed quantities 
with the anisotropic spacetime results in the ampli- 
fication of the convective field 5E\ ~ ^SB{t)5u'^ , 
which is the only one to survive in the ideal-plasma- 
limit [e.g. see Eq (30)]. Through Faraday's law, 
any amplification in the convective field leads to an 
analogous growth in SB'^ , at the expense of the cos- 
mological expansion. Accordingly, after exhausting 
any available amount of energy, the magnetic field 
perturbation reaches at a maximum value, before 
it is suppressed due to the cosmological redshift. 

• On the other hand, for 1] ^ 0, the magnetic field 
perturbation also increases rapidly at early times 
after Iq {At ~ 10^^ sec), reaching at values up to 3 
times its initial one. However, in this case, the per- 
turbation's amplitude is saturated, acquiring suffi- 
ciently large values for long enough time intervals 
[At ~ 10^^ sec). This is due to the fact that, be- 
sides the convective field SEf, a non-zero resistiv- 
ity favors also convective currents {SE2 = r/ SJ^) 
(For ?? ^ 0, the Ihs of Eq (30) corresponds, through 
Ampere's law, to an electric current). Accordingly, 
now, the available energy amount to be absorbed 
by the perturbed quantities is larger, and therefore, 
the magnetic perturbation remains at high levels 
for longer time intervals. 

As a result, after saturation, the magnitude of SB 
constitutes a fraction of 5 x 10^"^ of the unperturbed 
value of the magnetic strength. In this case, the 
quadrupole anisotropy induced in the MRB along 
the X-axis, reads 



AT 



rj2 



SB, 
B ' 



(77) 



resulting in 

AT 

— ~ 4.10 X 10"^ (78) 

X 

i.e. the corresponding value is enlarged by 1%. 



The numerical results indicate a completely different 
behavior for the electric field perturbation (Fig. 2). Not 
only the growth rate and the highest value of 5Ey are a 



03 
03 



to 




1.05 



t(x10^^sec) 



FIG. 1: The time-evolution of the magnetic field perturba- 
tion, for Bo = 7 gauss and for several values of the resistivity 
77 (sec). Notice that, for 7^ 0, the perturbation's ampli- 
tude is saturated, acquiring large values for long enough time 
intervals. 



little bit smaller than the corresponding values of (5_B^, 
but, also, the suppression rate of the perturbation's am- 
plitude, is much larger than that of 5B^ ^ resulting in a 
rapid decrease of the electric field at late times. It ap- 
pears that the expanding Universe disfavors strong elec- 
tric fields. 



n=0.0870 
n=0.0745 
n=0.0645 
n=0.0000 




o 
> 

2. 0.018 
llT 

5^ 0.016 



t(xlO^^sec) 



FIG. 2: The time-evolution of the electric field perturbation, 
for Bo — 7 gauss and for several values of the resistivity 77 
(sec). Notice that, in this case, there is no saturation. 



We conclude that, for reasonable values of the resistiv- 
ity, the magnetic field perturbations are lead to a real in- 
stability, acquiring large values for sufficiently long time- 
intervals. The influence of resistivity, in triggering insta- 
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bilities in anisotropic cosmological models, has been the 
subject of research in the past (Fennelly 1980). To the 
best of our knowledge, however, this is the first time that 
a direct connection between the resistivity and the satu- 
ration of the perturbations' amplitude at high values for 
long time intervals, is suggested and discussed. 



VII. JEANS-LIKE INSTABILITIES 

The question that arises now is, if the cosmological 
model under consideration admits also other kinds of in- 
stability and which is their role in connection to the resis- 
tive one. To answer this question, we study the evolution 
of purely gravitational perturbations, examining whether 
they admit a growing behavior (Jeans instability) or not. 

In the absence of e/m fields (and their fluctuations), 
one is left with the system of the perturbation equations 



{Hr + 2Hs - iLu)Sp = -ikp{t)Su^ 



(79) 



im + [Hr + 2Hs - iLu)p{t)]Su' = -i—C^Jf(80) 
the combination of which, yields 

[^^ + iHn + 2Hs)^iLu][iHR + 2Hs)-iu;] = -^C^ (81) 

Taking into account the particles' number conservation 
law, Eq (81) results in 



iu;{Hr + 2Hs) ~ kly^Cl 







(82) 



which, with the aid of the corresponding Fricdmann 
equation H"^ = ^^p, reads 



52 ■ 



Eq (88) is identical to the isotropic (FRW) result, pre- 
dicted by Weinberg (1972), in the relativistic theory of 
small fluctuations. 

In contrast to the high frequency e/m waves (52), as 
regards the corresponding energy-density perturbations, 
propagation is possible only when their physical wave- 
number is larger than a characteristic value, arising from 
the condition (84), otherwise, after some time they be- 
come unstable and grow exponentially with time (Jeans- 
like instabilities). 

Taking into account the background solution (15), Eq 
(84) at t = to reads 



1 



Csto 



1.86 X 10"^^ cm-i 



and the corresponding Jeans length is given by 
Ac = 27r Cs to = 3.37 x 10^^ cm 



(89) 



(90) 



Propagation of density perturbations with A > Ac is not 
possible, for every t > to and we are lead to a gravita- 
tional instability. The larger the coordinate wave-length 
is, the most prominent the unstable behavior will be (Fig. 
3). 



describing damped acoustic waves. With respect to lo, Eq 
(82) is a second order algebraic equation with roots 



.Hr + 2Hs 



2Hs 



provided that 



kphCs ^ 



Hr + 2Hs 



Y (83) 



(84) 



In this case, the energy-density perturbations (36) reduce 
to 



5p. 



A 



sVr 



where, cu^i is given by 

fc2 



H„ + 2H 



S\2 



(85) 



(86) 



Eq (86) represents the dispersion relation for the prop- 
agation of the energy-density fluctuations. It is worth 
noting that, in the isotropic case, where Hr = H = Hs, 
it yields 



LUr 



C 



9 



H' 



(87) 



o 

(D 



Q. 

to 




t(x10^^sec) 



FIG. 3: The time-evolution of the energy density perturba- 
tion, for several values of the coordinate wave-length A in 
terms of Ac. 

Furthermore, one may ask whether the waves with 
wave-number around kc in a non-ideal plasma, may grow 
faster than those in an ideal plasma. For every t > to, 
the physical Jeans length along the i-axis [A^; — XcR(t)], 
is larger than the corresponding length along the other 
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two axes [Xy = XcS{t) = X^], due to the background 
anisotropy, suggesting formation of "cigar-like" conden- 
sations within the anisotropic fluid. Since this fluid is 
conductive, these condensations act in favor of electric 
currents which may lead to a further amplification of the 
e/m perturbations, fortifying any pre-existing resistive 
instability (Fig. 4). Therefore, a Jeans-like instability 
enhances the phenomena related to the resistivity. 



0.035 



03 
03 



to 




1.01 



1.02 



1.03 



1.04 



1.05 



t(x10^^sec) 



FIG. 4: The time-evolution of the magnetic field perturbation 
for rj = 0.0645 sec, Bo = 7 gauss and for several values of the 
coordinate wave- number in terms of kc- 

On the other hand, numerical results indicate that 
waves with wavelength around kc become more promi- 
nent as the resistivity grows (Fig. 5). This result also 
has a clear physical interpretation: As we have already 
seen, any increase in the resistivity fortifies the surround- 
ing magnetic field. A strong magnetic field organizes 
plasma along its lines, favoring any pre-existing conden- 
sations. Hence, resistive instabilities act in favor of the 
corresponding gravitational ones and vice versa. 



VIII. DISCUSSION 

We study the evolution of the magnetosonic waves in a 
magnetized, resistive plasma, which governs the dynam- 
ics of an anisotropic cosmological model. After construct- 
ing the general set of MHD and Einstein's equations for 
the anisotropic cosmological model (see the Appendices), 
we solve the field equations to obtain the zcroth-order 
solution. In order to determine the waveforms admitted 
by this system in the first place, we introduce wave-like 
perturbations and, neglecting all terms higher or equal 
than the second order, we extract the dispersion rela- 
tion at t = to, i.e. at the beginning of the interaction 
between magnetized plasma and curved spacetime. It 
appears that magnetosonic modes can be excited due to 
the anisotropy and the resistivity. For t > Iq, we in- 



o 

CD 



Q. 

!0 



0.01415 

0.0141 
0.01405 
0.014 
0.01395 

0.0139 
0.01385 

0.0138 
0.01375 

0.0137 
0.01365 



n=0.0870 
n=0.0745 
n=0.0645 



1 1.02 1.04 1.06 1.08 
t(x 10^^ sec) 



1.1 



FIG. 5: The time-evolution of the energy density perturba- 
tion, for k = kc and for several values of the resistivity rj 
(sec). We observe that the Jeans instability becomes more 
prominent as the resistivity grows. 



tegrate numerically the perturbed equations, using the 
dispersion relation for the fast-magnetosonic waves. 

We find that, at early times, both the electric and the 
magnetic field perturbations grow exponentially, at least, 
in the regime where the linear analysis holds. However 
there is a major difference in their behavior in the pres- 
ence of a non-zero resistivity. For 77 0, the magnetic 
field perturbation after increasing to reach at values up 
to 3 times its initial one, is subsequently saturated, re- 
maining at high levels for sufficiently long time intervals 
{At - 10^2 sec). 

The situation is completely different as regards the 
electric field perturbation. Not only the growth rate 
and the highest value of SE^ are a little bit smaller than 
the corresponding values of but, also, the suppres- 
sion rate of the perturbation's amplitude, is much larger 
than that of SB^ . Accordingly, the electric field decreases 
rapidly at late times. It appears that the expanding Uni- 
verse disfavors strong electric fields. 

Finally, we have shown that waves with wave-number 
around kc, are enhanced in non ideal plasmas. 
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Appendix A 



We present the closed set of MHD and Einstein equa- 
tions (in the system of units where ft = 1 = c) for the 
anisotropic cosmological models of Bianchi-Type I 

ds^ = -dt^+R^{t,z)dx^+S^{t,z)dy^+T^{t,z)dz^ (Al) 

in the presence of an anisotropic perfect fluid, which al- 
lows for acoustic waves along the z-direction 



Pig 



and an e/m field of the form 



{) {) E y 



-Ey 

-B^ 



(A3) 



where, Greek indices refer to the four-dimensional space- 
time and Latin indices refer to the three-dimensional spa- 
tial section. In Eq (A3), F"'^ is the Faraday tensor in flat 
spacetime. The components of the e/m field in curved 
spacetime are defined by 



(A4) 



where, the non-zero components of the orthonormal 
tetrad of the local Lorentz frame for the metric (Al), 
are given by 

ef = (-1, 0, 0, 0) el = (0, 0, 0) 



1 



R S 



ST T R 



T2 \ RS 



ST 



TR 



^ ^TiGp + G[{Eyf + [B'^f] 
The (xx)-component is given by 



(A7) 



T 



i?2 S 



ST\ R^ T S 



T2 \ S 



ST 



T2 \ TS 



-S 



= SirGR^P,, - GR''[-{Eyf + [B^f] (AS) 
The (yy)-component is given by 

'r f\ S'^fR"\ „^ f Rf\ S^(t'r^\ 



R T T^ \ R 



-S^ 



RT 
RT 



T2 \ TR 



= 8TTGS^Py + GS'^[-{Ey'f + {3"=^] {A9) 
The (zz)-component is given by 

\ S R ^ T ' \RS \ RS 



e^- (0,0,1,0) = (0,0,0,1) (A5) 

Therefore, in the curved spacetime (Al), the non-zero 
components of the Faraday tensor are 



^02 ^ 



-F 



20 



= SnGT'^p, + GT^[+{Eyf + [B^'f] 
The (tz)-component is given by 

RT 



R 

"r 



5^ 



RT 



S T 
'ST 



-SttGT^ 



pu 



iGTEyB"" 



(AlO) 



{All) 



77123 

^ ST^ 



In what follows, the dot denotes time-derivative, while 
the prime denotes differentiation with respect to z. The 
Einstein equations, G^^ = 87rG(7^{f"''' - 
The (tt)-component is given by 



Appendix B 

On the other hand, the Maxwell equations in curved 
spacetime are written in the form 




ST TR 
'sf^TR 
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where, J" = ^F^^Ufs is the current density and r] is the 
electric resistivity of the fluid. Accordingly, we obtain 



-dtsy 




R' 



and 



S' 



~d,{SEy) + dt{STB'') =0 



Appendix C 

Taking the time and space component of = 0, we 
obtain the required equations of motion in a covariant 
form, namely 



1 



1 



dtW) + ^{E^ + B^)] + d^pu' - —EB 



Stt 



47rT 



'T3 



p,{t) + -{E^+B^) 




pu 



AttT 



-EB 



R 
R 



Piit)--i-E' + B^] 



P2{t) 



Stt 



{-E^ + B^) 





S' 




i- 








2^2 



T 
T 



psit) + —i+E' + B^] 
an 



P3{t) + —{E^ + B^) =0 

In addition, the particles' number conservation law 
{pu^);n = 0, results in 

dtipu'^) + dAfiu')+ 




p{t) + -{E^ + B'~) 
an 



(CI) 



^EB]+dA^[p,it) + l-iE^ + B^)]y 




(C3) 
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